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We study the applicability of magnetization and specific heat equations derived from a lowest- 
Landau-level (LLL) calculation, to the high-temperature superconducting (HTSC) materials of the 
YBa2Cu307_5 (YBCO) family. We find that significant information about these materials can 
be obtained from this analysis, even though the three-dimensional LLL functions are not quite as 
successful in describing them as the corresponding two-dimensional functions are in describing data 
for the more anisotropic HTSC Bi- and Tl-based materials. The results discussed include scaling 
fits, an alternative explanation for data claimed as evidence for a second order fiux lattice melting 
transition, and reasons why 3DXY scaling may have less significance than previously believed. We 
also demonstrate how 3DXY scaling does not describe the specific heat data of YBCO samples in 
the critical region. Throughout the paper, the importance of checking the actual scaling functions, 
not merely scaling behavior, is stressed. 



I. INTRODUCTION 

The critical behavior of the high-temperature super- 
conductors in finite magnetic fields applied perpendicular 
to the copper oxide planes has been described by both 
lowest Landau level (LLL) theory and three-dimensional 
(3D) XY theory with varying degrees of success. It is 
widely expected that 3DXY behavior should hold at low 
fields and that LLL should be valid at higher fields. There 
is however little consensus about what the value of the 
crossover field should be and how well either of these 
theories describe resistivity data, magnetization data, or 
specific heat data pj-H]. One group |p[ has claimed that 
LLL should not be valid for fields less than ten Tesla 
(T) in deoxygenated YBa2Cu307_5 (YBCO) thin films 
based on conductivity measurements while three of the 
present authors P,R| have found it to be valid down to 
approximately two Tesla based on an analysis of specific 
heat data from YBCO and LuBa2Cu307_5 (LBCO) sin- 
gle crystals. 

On the theoretical front, considerable work has been 
done-namely, the derivation of analytic expressions for 
the magnetization and specific heat LLL scaling functions 
for two dimensional (2D), three-dimensional (3D), and 
layered systems. This was achieved [g|,0 by using the LLL 
approximation in the Ginzburg-Landau (GL) formalism. 
The two-dimensional portion of this work has had strik- 
ing success in describing the magnetization of the highly 
anisotropic high-temperature superconducting materials 
BiaSraCaCusOg (BSCCO-2212) and the Tl-based com- 
pounds for magnetic fields applied perpendicular to the 
copper oxide planes. For example, the two-dimensional 
(2D) functions HJ^ have a field independent value at a 
particular temperature T* pC| ] which is also a crossing 
point for the magnetization curves. Such behavior has 
been observed in BSCCO-2212 by many authors (ll| and 



furthermore, Wahl and coworkers |12|] have not only ob- 
served such crossover in their magnetization data from 
Tl-based single crystals but they have also fit the 2D 
functions of Ref. |g,|[ to their data finding good agree- 
ment. 

Little work has been done to fit the theoretical func- 
tions to specific heat data on the highly anisotropic 
HTSC materials. Kobayashi et al . |13[ are among the 
few to publish specific heat data |14[| for various fields 
near the critical temperature on such compounds. They 
scaled their specific heat data from a c-axis aligned 
(Bi,Pb)2Sr2Ca2Cu302; bulk sample and compared it to 
the 2D scaling function of Tcsanovic and coworkers [pP] 
finding reasonable agreement. They also found a crossing 
point in their magnetization data. 

Even less has been done to compare the theoretical 
expressions and scaling functions to experimental data 
for the more isotropic YBCO materials. In Ref. |^, an 
approximation to the 3D LLL specific heat function was 
compared to scaled specific heat data from various YBCO 
samples (including a YBCO single crystal from Ref. Q|) 
and a LBCO sample, with satisfactory agreement. Fur- 
ther, in the work of Ref. [J] , although a quantitative com- 
parison was not made, one can find qualitative agreement 
between the scaled temperature derivatives of the specific 
heat and the second derivative of the 3D magnetization 
function. (See Eq. (5) and Figs. 6 and 8 of that reference.) 
Lastly, we are not aware of any work comparing the 3D 
theoretical magnetization function to magnetization data 
of YBCO-class materials. This can be attributed in part 
to the complexity ||l^ of the 3D specific heat and mag- 
netization functions of Ref. l9l . 

In this paper, we examine the 3D specific heat and 
magnetization functions of Ref. 0] comparing them to 
data from YBCO samples. Not only is such a compari- 
son lacking and certainly needed in order to learn more 



about the validity of the theory but we will see that it 
yields valuable insights into other questions about the be- 
havior of the YBCO materials, besides the nature of their 
fluctuation behavior. For example, we will demonstrate 
how data presented as evidence for a second order flux 
lattice melting transition can be explained within LLL 
theory without invoking flux lattice melting arguments. 
Furthermore, it will also be seen that the 3DXY scal- 
ing is so general as to describe "theoretical data" derived 
using the LLL theory. This exemplifies the importance 
of knowing the actual scaling function and leads one to 
question the significance of 3DXY scaling. The failure of 
3DXY scaling to describe specific heat data from YBCO 
samples will also be shown. 

Our focus here will not be so much on the validity 
of LLL scaling for the HTSC materials as on the appli- 
cability of the specific expressions and scaling functions 
from one particular calculation based on a nonpertur- 
bative approach to the GL-LLL theory |gj^. There are 
two separate issues here. First, there is the question of 
describing the HTSC's in a GL-LLL formalism, which 
has already been answered, in our opinion, through the 
success of LLL scaling. Second, while the expressions we 
use are of compact form and should be useful for analysis 
of experiments and phenomenology, they are not mathe- 
matically identical to the exact solution of the GL-LLL 
theory. Therefore, there is a need to address any possi- 
ble disagreement between theory and experiment arising 
from the additional approximations involved in the ex- 
pressions for the scaling functions of Refs. [pP] relative 
to the exact answer within the GL-LLL theory. Although 
this issue is essentially resolved for very anisotropic, "al- 
most" 2D HTSC systems where the 2D LLL scaling func- 
tions of Refs. BM are known to be very accurate, |16| it 
has not been investigated for the relatively isotropic ma- 
terials of the YBCO class. 

In this work, we take advantage of the availability of 
numerical work on the (quasi) 3D GL-LLL model [[l7[ . 
The existence of this numerical work for the magnetiza- 
tion will allow us, as we shall see, to determine some of 
the fitting parameters in a way that it is not constrained 
with experimental uncertainties involving, for example, 
the subtraction of "background" terms. We can say that 
we use these numerical results to "calibrate" certain pa- 
rameters in the scaling functions. This is very convenient, 
since the increased complexity of the 3D functions, as 
opposed to the 2D case, would otherwise make our task 
much more intricate and the conclusions weaker. 

The paper is organized as follows: The theoretical 
functions calculated from the nonperturbative approach 
Mm to the GL-LLL theory will be set forth and dis- 
cussed in Section II A . In Section LIB, the "calibration" 



fits of the numerical 3D magnetization data to the the- 
oretical result, [Eq. p)] are performed, and then fits to 
actual magnetization data from YBCO samples are done. 
Then, in pQ, we give a simple explanation of the pecu- 



liar behavior of the field dependence of the partial deriva- 
tive dM{H, T)/dT found |l§ in YBCO and BSCCO. We 
show that this behavior is simply explained in terms of 
the LLL scaling functions. Fits of the theoretical 3D 
specific heat function [Eq. (p[)] to specific heat data from 
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the same materials as in II B| are reported in Section II D 



along with an alternative explanation to data claimed as 
evidence for a second order flux lattice melting transition. 
Finally, implications of this work for 3DXY scaling and 
the importance of the scaling functions is demonstrated 



in Section [I E followed by a discussion and summary in 
Section f!I 



II. THEORETICAL FUNCTIONS AND DATA 
FITS 

A. Theoretical GL-LLL Functions 

As mentioned above, the 2D functions of Refs. 
have had considerable success in describing magnetiza- 
tion data from the highly anisotropic HTSC materials. 
Here, we will focus on the 3D specific heat and mag- 
netization functions since the 2D functions are already 
examined and the quasi-2D functions are less tractable. 
The magnetization as a function of applied magnetic field 
H and temperature T is written as Eq. (26) of Ref. ||] 
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The specific heat function can be rewritten (including 
important subleading terms) by considering |19[ the pure 
3D limit of the quasi 3D result (Eq. (30) of Ref. [||): 
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Recall that g is related to the temperature through 
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In the above equations, |^^ H'^2 is the first derivative of 
the critical field Hc2{T) with respect to temperature (and 
is assumed to be a constant), k — Xab/£,ab (where Aob and 
^ab are the penetration depth and coherence length re- 
spectively in the ab plane), ^ is ^c to within a multiplica- 
tive constant (^c is the coherence length along the c-axis), 
Tc{H) is the finite field mean- field transition tempera- 
ture, Tco — Tc{Q), Cmf{T) is the mean-field specific heat, 
4>o the superconducting flux quantum, ks the Boltzman 
constant, and Q, K^ and M are adjustable parameters 
whose values are roughly tt, \/2 and 2-\/2 respectively. 
The factor of 47r in Eq. (Q) is needed to convert M{H, T) 
which is in units of emu/cm"^ to Gauss (G). 

The fits to Eqs. (|l|) and (0) are nontrivial since U{G) 
[Eq. (||)] and G{g) [Eq. (^] are functions of one another 
and the temperature T is related to g through a transcen- 
dental equation. In place of a commercial fitting package, 
we have written a code which uses IMSL routines to fit 
the functions [Eqs. (|l|) and @\ (which must be calculated 
self-consistently) to the data. 

The 3D scaling functions [Eqs. (|l|) and (||)] are more 
complex than those of the 2D functions because there is 
an extra length scale which describes the bending of the 
vortex lines along the field direction (here taken to be the 
c-axis). For this reason, as we noted above, we will take 
advantage of the numerical results [|^ on the GL-LLL 
theory for the magnetization, which is the only quantity 
presently available from numerical work, to "calibrate" 
certain parameters in our 3D scaling functions. 



B. Magnetization 

As explained above, we will begin by narrowing down 
the number of parameters available to perform fits to ac- 
tual experimental data by first considering fits to results 
from a numerical calculation which simulates YBCO. We 
will then use our results from this fitting as a means of 
"calibrating" Eqs. (|I]) and §). 
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FIG. 1. (a) A fit of Eq. ([Z|j to the "numerical" magnetiza- 
tion data of Ref. llljj . (h) The same data |j'/|/ scaled according 
to 3D LLL theory. 

The numerical calculation to which we refer above was 
done by Sasik and Stroud [jl^l using a GL-LLL formu- 
lation for a layered system (with parameters similar to 
those of YBCO H'^^ = l-S?", « = 52, Tco = 93, and an 
anisotropy factor 7 — ^ab/S.c = 5) in order to study flux 
lattice melting. Their results are reproduced in Fig. ||a. 
We have done a four-parameter fit of Eq. (|l|) to this nu- 
merical "data" |l^] and plot our results as lines (solid for 
2T, dashed for 3T, and dotted for 5T) with the data in 
Fig. Ha. As one can see, the agreement is excellent. The 
fitting parameters are Q, K, M, and the constant A relat- 
ing S, to Cc and we find Q = 10.25, K = -5.95, M = 7.38, 
and ^ = 0.2918^c (i-e., A = 0.2918) since £,c = 2.821 here. 
To verify the consistency of the numerical data with LLL 
theory, we scaled the data according to the 3D LLL form 
M{H,T)/{HTf/^ = f{[T-T^{H)]l{THf/^) and found 
the data to collapse flawlessly as shown in Fig. |l|b. When 
doing such scaling, one can typically use Tc{H) as an 
adjustable parameter (and, to a lesser extent, the back- 
ground parameters), but since H'^2 and Tco are known 
in this case, there are no adjustable parameters, which 
makes the scaling of this data very convincing. 

Using these values of Q, K, M, and A determined 



from the fit to the numerical "data" we then performed 
a fit to the two and three Tesla magnetization data 
of Jeandupeux et al |^1|]. There are nine fitting pa- 
rameters, namely H'^^, '«: Cc Tc{2T), Tc{3T), Bo{2T), 
Bi(2T), -Bo(3T), and Bi(3T) where BoiH) and Bi{H) 
are the field-dependent constants used to adjust the sub- 
tracted background: Mb = iBo + Bi/T)H |22). Three- 
parameter fits were then used to find Tc{H), Bo{H), 
and Bi{H) for the four and five Tesla fields. We find 
fc = 3.781, H'^2 = 1-S37T/K, k = 56.02, and TdH) = 
90.917^, 90.38K, 89MK, and 89.04is: for H = 2T, 3T, 
4:T, and 5T respectively and show the fits in Fig. 0. The 
fit to the 5T data is the least satisfactory which we at- 
tribute to the data: One can see from Fig. 3 of Ref. ||2^] 
that the 4T and 5T data have spurious behavior at low 
temperatures instead of collapsing to the mean-field tem- 
perature dependence which could be a result of the entry 
into the irreversible region. Except for the 5T fit, the fits 
are reasonable and the parameter values are similar to 
those found by others which gives us confidence that the 
theory is a good description of the data. Furthermore, 
when one compares the value of H'^2 obtained in the fits 
to the values H'^^ = 1.8?>T/K found from the Tc(i7)'s 
for H = 2-4T (throwing out the less satisfactory 5T fit), 
one finds good agreement strengthening the credibility 
of the fit. If the quantities Q, K, and M are added as 
fitting parameters, rather than being taken as obtained 
from the numerical "calibration" , significantly better fits 
are not obtained. This is as expected, if the procedure is 
correct. 

In the Inset to Fig. 0, we show the LLL scaling of the 



magnetization data of Ref. |21| using the parameters ob- 
tained from our fit. The collapse of the data is good 
in the critical region but fans out somewhat in the low- 
temperature, mean field region. We believe that some 
of this fanning may be due to the spurious behavior as- 
sociated with the irreversible region which we discussed 
above. [g3[ 

We have also attempted fits to the magnetization data 
on a YBCO single crystal by Salem-Sugui and da Silva 
^4| . In this case the fits were rather poor unless unphys- 
ical parameter values were chosen. The situation did not 
improve when Q, K, and M were added as fitting pa- 
rameters. This is in contradistinction with what occurs 
with the numerical results and with the data of Ref. [gl| . 
We are more inclined to believe that the origin of the 
discrepancy lies with the data set than with the theory, 
but we cannot be certain until more magnetization data 
from YBCO becomes available to us. 

We now turn to using Eq. p) to explain the results 
of Ref. |l8[ in which it was found that the temperature 
derivative of the magnetization has an approximately 
field independent value and a crossover point at a dis- 
tinct temperature whose value lies close to that of Tco 
for YBCO. 
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FIG. 2. The magnetization data from a YBCO twinned 
single crystal published in Ref. fej/ along with the fits to 
Eq. in). [INSET: Our LLL scaling of the magnetization data 
of Ref. \2aI shown over a larger range than that of the authors 
of Ref. Wnj which is in the Inset to Fig. 5 of that reference. 
The y-axis is A-KM/{THf^^ x lQ'^[G/{OeKf'^ and the x-axis 
IS t X lQ^[K^'^/Oe^'^] where t is defined m Eq. |./ 



Behavior of the Magnetization Temperature 
Derivative 



The behavior of the temperature derivative of the mag- 
netization dM{H,T)/dT has been studied as a function 
of the field H in Ref. |l8| . The behavior is quite strik- 
ing. The most salient feature of the experimental results 
(Fig. 2 of Ref. |18|) is the very weak dependence on the 
field of this partial derivative at temperatures near the 
mean-field temperature. This weak dependence extends 
to a rather wide temperature range (more than eight 
degrees) and in fact, can nearly be called a field inde- 
pendence for the fields 7J > 2 for both the YBCO and 
BSCCO data. Also remarkable is the feature upon which 
the authors in the cited experimental work focused which 
is the apparent crossing of the data at a temperature very 
close to Tco- 

These experimental results can be easily understood 



from the LLL scaling formula, Eq. (|l|). This equation 
can be rewritten simply as: 



M = {HTf/^iL{t) 



(8) 



where /x is the scaling function and t the scaled temper- 
ature variable defined in Eq. (R). From these equations 
we have for the temperature derivative: 

where ^'{t) is the derivative of /i(i) with respect to its 
argument. We can the proceed to the evaluation of the 
mixed second partial derivative, with the result: 

rP'Af A 1 
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where we have used dTc{H)/dH = l/H'^2- Because the 
scaling function for the magnetization is essentially a lin- 
ear function of its argument (except in a very narrow 
region near its kink) , we can drop the term proportional 
to /i" in Eq. (|lO|). We then have that the condition for 
the mixed derivative to vanish is. 



(i?T)"V3^(t) = M'(t)^^^|^, 
which, taking into account Eq. (|7|), can be written 
^j.{t) = ^'{t)t w ^'(0)i. 



(11) 



(12) 



Since /l((0) is small, this relation is to good accuracy the 
Taylor expansion of /i(t) about t = and it will triv- 
ially hold over an extended region. Indeed, since fj, is 
nearly everywhere linear, so that there are no higher or- 
der terms, it may appear that we have nearly proved that 
the vanishing of the mixed derivative is an identity. This 
is not at all the case, chiefly because /i(0) cannot ev- 
erywhere be neglected. However, the argument makes it 
abundantly clear that the weak dependence on the field 
of the mixed derivative of AI and the crossing point follow 
easily from LLL scaling. 



D. Specific Heat 

In this section we will examine the specific heat func- 
tion [Eq. (|)] fitting the data of Ref. ||l| to it. We will 
then use Eq. (0) to address the results of Ref. ES] in 
which evidence for a second order flux lattice melting 
transition is claimed and show how such experimental 
results can be alternatively explained via Eq. (0) with- 
out need to invoke flux lattice melting arguments. 
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FIG. 3. Theoretical specific heat curves for various fields 
(H — l-ST plotted at 0.5 T intervals) calculated from Eq. M). 

We begin by examining the specific heat curves pro- 
duced by Eq. (g). We display this function for various 
fields in the vicinity of the critical area using parame- 
ters characteristic of YBCO in Fig. ^. For Cmf{T), we 
used a standard form: Cmf{T) = iT[l + b{T/Tco - 1)] 
which was taken to be a constant above Tc{H). This 
form produces an artificial nonanalyticity of the curves 
at Tc{H) since the first derivatives are not continuous 
but this is unimportant since it does not affect our anal- 
ysis in any significant way. As one can see in Fig. 0, the 
curves generated from Eq. (0), which have the features 
of a mean-field "ramp" with fluctuations which produce 
a peak, are qualitatively similar to the data from YBCO 
samples. See, for example. Fig. 2 of Ref. |2^] (reproduced 
here in Fig. |), Fig. 1 of Ref. [||, or Fig. 2 of Ref. 0. 
What is common to all of these curves is that they do not 
collapse immediately for temperatures below the peak. 
This is in agreement with a result derived from LLL the- 
ory, as originally pointed out in Ref. [|7|- There, using 
the Maxwell relation, {d^M/dT^)H = {dCH/dH)T/T, it 
was noted that the left hand side of this equation is pos- 
itive for lower temperatures, which means that C{T, H) 
must increase with increasing field |26|. As we will dis- 
cuss below (Section HE), 3DXY theory cannot account 
for such behavior. 

We have attempted fits of the specific heat data of 
Ref. Q to Eq. (|). The best test of the theory would be 
to be able to fit the specific heat data with the param- 
eters obtained for the magnetization data on the same 



sample, which was discussed in Section HB. This could 



not be done however, and a good fit to the entire rele- 
vant temperature range could not be obtained even when 
lifting the constraints from the magnetization fits. We 
believe this is due to the 3D specific heat function above 
the peaks not having the same qualitative behavior as 
data from YBCO samples which we believe exhibits 2D 
behavior in this region. The experimental side of this 
statement is reasonable since 2D fluctuations have been 



observed in these materials through electronic transport 
measurements |g^. (For evidence in specific heat data, 
see Refs. |g,0-) The theoretical side of this statement is 
plausible for the following reasons. As one can see, the 
theoretical curves for three dimensions clearly exhibit a 
crossover point above the peak. While this crossover can 
be shifted to higher temperatures by varying the parame- 
ters, it still tends to occur at smaller temperatures than it 
does for the 2D theoretical curves. More important, the 
2D curves decay more quickly to zero above the transi- 
tion. 
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FIG. 4. YBCO specific heat data for various fields 

(H = 3-7T) from Ref. /||/ along with the fits to Eq. (gj. 

To account for the apparent dimensional crossover of 
the specific heat data, we have done fits of the 2D func- 
tion [g] to the data omitting a sizeable temperature win- 
dow [83.3K:91.5K] around the peak where the 3D behav- 
ior is expected to dominate in order to fix the background 
{Cb = (-8.6336 + 2.2602T) mJ/gK). We then do a fit of 
the 3D function to the 3T, 4T, 5T, 6T, and 7T specific 
heat data for the temperature range [82K:89K]. The re- 
sults are good as can be seen from the parameter values: 
H'^2 = 1-82T/X, K = 65.6, and ^^ = 3.1819. This is re- 
inforced by a visual inspection of the fit in Fig. where 
the curves agree quite well in the region we believe to be 
3D. One also notices that the curves are on the high side 
at temperatures above the peak, which agrees with our 
earlier statement that the 2D fluctuations decay more 
quickly with temperature than do the 3D ones. 

The parameter space is large here and it would be im- 
possible to explore all of it to determine the very best 
fit. From our investigations however, we are certain that 
an improved fit with better parameter values can be ob- 
tained for example by allowing Q, K , and M to vary. One 
could also allow for a quadratic term in the background or 
even in the mean-field term which we have extended over 
a large temperature range. Even without exploring the 
large parameter space, we have demonstrated the agree- 
ment between the 3D LLL specific heat function [Eq. (||)] 



and YBCO data. 

In the remainder of this subsection, we will discuss 
more indirect consequences of LLL theory as applied to 
the specific heat of YBCO class materials. We will show 
how using Eq. (0) we can give an alternative explanation 
for features that have been claimed p5[ as conclusive 
evidence for second order flux lattice melting without in- 
voking flux lattice melting arguments. As we pointed out 
above, the curves in Fig. produced from Eq. (0) repro- 
duce the key features of specific heat data from YBCO 
samples, such as that seen in Fig. 1 of Ref. ||2^], especially 
for temperatures below the peak. As the field increases, 
the peak in the theoretical curves moves down in tem- 
perature and broadens. It is this general behavior which 
we claim produces the feature that Roulin et al. cite as 
evidence for second order fiux lattice melting. 

To arrive at their conclusion of second-order flux lat- 
tice melting those authors [^ subtract a data set for 
one field from a data set from a slightly larger field and 
find a peak in the "differential" . This peak is due to the 
higher field starting to peak at a slightly lower tempera- 
ture due to peak broadening and transition temperature 
suppression. To see if this peak corresponds to a flux 
lattice melting, they look at the field dependence of the 
peak temperature and find that it agrees with the same 
curve found in Ref. p8| by analysis of magnetization and 
resistivity curves to derive a flux lattice melting line for 
a YBCO sample. Although this seems persuasive, we 
will now show that similar results can be obtained from 
Eq. d). 

Consider Fig. |a. Here, SC{T,H = 4.25T) = 
C{T, H = 4.5T) - C{T, H = 4.0T), the difference of two 
theoretical curves calculated from Eq. (0) , is plotted and 
is seen to have the same behavior as that of Fig. 2b in 
Ref. [^: a peak followed by a deep trough. The field 
dependence of the temperature at which each SC{T,H) 
peaks is then calculated. This is plotted in Fig. ||b along 
with two lines. The first line (solid line) is the best three 
parameter fit of this theoretical LLL result to the func- 
tion H{T) = a{l ~ T/hf with a = 111.2, b = 92.1 
and c — 1.48 (standard deviation= 0.02). The second 
(dashed) line is a two parameter fit with a fixed c — 1.33 
and a = 89.2 and h = 91.1 (standard deviation= 0.05). 
One can see that the while the exponents differ by 11%, 
the two curves are hard to distinguish. These lines are 
to be contrasted to that of Refs. p5| , p8| where the corre- 
sponding line, H^{T) = 99.7[l-r/92.5]i-3^ is identified 
with the second order melting line. The numbers in both 
fits agree quite well with each other, especially the expo- 
nent which is the most critical parameter in this curve. 
However, in our case it obviously has nothing to do with 
mehing. The curve H„(T) = 99.7[1 - T/Tc[0)Y-^^ has 
not since been experimentally reproduced with the same 
numbers. It therefore seems that the agreement between 
the specific heat data of Ref. pS] and that of Ref. pa] 
may be just coincidental, and reflect a property of the 



LLL specific heat. 
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FIG. 5. (a) A "differential" (SC{T,H = 4.25T) = 
C{T, H = 4.5T) - C(T, ii' = A.OT)) from two of the theo- 
retical curves in Fig. ^. This curve has the same qualitative 
behavior as Fig. 2b of Ref. feq/. (b) The field value versus the 
temperature at which the "differentials" peak (as explained in 
the text) for those field s and two theoretical fits which are dis- 
cussed in Section HE . 



The authors of Ref. |2^] have since come out with sep- 
arate and more convincing evidence p9| for a first order 
phase transition in a cleaner YBCO sample finding data 
similar to that of Schilling et al |30|. Fig. 2b of Ref. [|9), 
the analog of their Fig. 2a in Ref. ]25|], shows clear spikes 
that were not present in their original paper. Our argu- 
ment above does not, of course, affect this independent 
evidence for a first order flux lattice melting transition 
which now seems to be well documented, f^,^ 



E. Implications for 3DXY Scaling 

In this subsection, we will address an important ques- 
tion concerning the experimental discrimination between 
LLL and 3DXY scaling. We will show that the LLL the- 
oretical results derived from Eqs. (|l]) and (g) can be used 



to generate curves which can then apparently be scaled 
in accordance with 3DXY theory. With such a result, 
one has to question the meaning of 3DXY scaling in ex- 
perimental data: finding 3DXY scaling does not at all 
exclude that the data actually is in agreement with GL- 
LLL theory. Further, the importance of comparing the 
experimental results with the actual functions is made 
obvious. This is a problem for 3DXY theory since no 
scaling functions are available for it. 

There are two scaling forms that have been used to an- 
alyze specific heat data according to 3DXY scaling the- 
ory. The first is that derived in Ref. p^: 



^0.0097 J^ ^0.747 



(13) 



where T^ is the zero-field critical temperature and the 
second is that derived in Ref. ||l|: 






^0.0097 



^0.747 



(14) 



where Co is the height of the specific heat cusp and the 
subscript SC signifies that it is only the superconduct- 
ing contribution, with the background subtracted out. 
(The exponent for H on the right hand side (RHS) of 
these equations is derived from the specific heat expo- 
nent a and we have used the value derived from He 
experiments. The theoretical value of a is 0.005.) The 
second scaling form m is the more convincing one since 
the scaling takes place over a wider and less trivial range. 
In the first form on the other hand, one is scaling mostly 
horizontal lines which are equal to zero [ p2[ , and thus 
are virtually guaranteed to scale. The first form thereby 
acts as a sufficient condition for 3DXY scaling and we 
use it here. If it does not scale according to this form, it 
certainly won't scale according to Eq. (Of). 

Thus, we proceed to generate "data" from the theoreti- 
cal expression Eq. (^). This is a portion of the theoretical 
results shown earlier in Fig. y. Then, we attempt to scale 
this "data" according to the 3DXY formula. The scaling 
results are shown in Fig. | for the 4T, 5T, 6T, 7T, and 
8T fields and as one can see, the collapse is reasonable. 
We have left out the smaller fields which as one would 
expect do not collapse onto these curves. If one were 
to consider the noise, and the background subtractions, 
which inevitably enters enter into the analysis of actual 
experimental data, one could call our scaling convincing. 
(When the exponent 0.005 is used for the exponent for 
H on the RHS of Eq. (|l|), the "data" does not collapse 
quite as well.) That LLL calculated "data" can be made 
to scale according to 3DXY theory exemplifies the exces- 
sive generality of 3DXY scaling and again raises questions 
on the significance of it. 
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FIG. 6. The theoretical specific heat curves from an LLL 
calculated Eq. m) (Fig. LSj scaled according to 3DXY theory: 
Eq.^. 

As a further insight on the problems with 3DXY scal- 
ing, we comment on actual specific heat data from a 
YBCO sample which was scaled using Eq. O the same 
way as the theoretical curves in Fig. pi This is the spe- 
cific heat data of Ref. |21^1, the same data to which we 



attempted fits in Section II D. The 3DXY scaling results 
in that work are shown in Fig. 4 in Rcf. [El|. The large 
central peak is the zero-field peak and so of course it is 
unimportant since it is simply an artifact of having sub- 
tracted off the zero-field data. The important region to 
consider is just to the left of the central peak, since that 
area represents the peaks (or the region of critical be- 
havior) of the finite-field data sets. As one can see, the 
data is not close to collapsing here. Because this data is 
plotted on a scale which brings in the size of the irrele- 
vant, large zero- field peak, scaling does not appear to be 
a dramatic failure. However when compared to the col- 
lapse of the data at temperatures above the peak where 
there are no fluctuations and one is dealing with only the 
background, one can see that the failure of 3DXY scaling 
is more obvious. We have done our own similar analysis 
on this data which we show in Fig. R in a more restricted 
scale. (We emphasize the point that if the data docs not 
scale with Eq. O, it will not scale with the more con- 
vincing Eq. n4[) It is seen that the 2T and 3T data do 
collapse but the IT data rides high and the larger fields 
go low. The failure of the scaled data to collapse in the 
actual critical region makes one question whether 3DXY 
behavior is valid. This is, again, one of the reasons why 
derivations of the actual functions like those of Ref. 0] 
are so important. 

We close this section on implications for 3DXY scaling 
with a discussion of the references to 3DXY in the flux 
lattice melting literature. There a melting line of the 
form H.JT)^ 99.7[I - T/Tc{Q)Y-'^'^ is found by many 
authors |^,^,^ who point out that the exponent here 
is the nearly the same as that (1.33) expected for 3DXY 



critical point analysis. And while some of the same au- 
thors [^ note that such an analysis appears incompat- 
ible with a first-order transition, that we get the same 
exponent from a LLL approach appears to make such an 
identification with 3DXY theory even more unlikely. 
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FIG. 7. The data of Ref. ^ scaled according to 3DXY 
theory. The critical region is to the left of the large peak which 
is the zero-field peak. One can see that 3DXY scaling does not 
describe the critical region. There are seven data sets plotted, 
corresponding to fields of 1 to 7 T at IT intervals. Higher 
field sets have lowers values at the extreme left of the graph. 



III. DISCUSSION AND SUMMARY 

In this paper, we have studied the relevance of 3D 
GL-LLL theory as exemplified by Eqs. dJ) and (H) for 
the magnetization and specific heat, to the relatively 
isotropic HTSC materials of the YBCO family. 

Eq. (w) was found to accurately describe the magneti- 
zation data for a YBCO sample |2^] and the numerical 
data modeling these materials in the GL-LLL formalism 
p7| . In the former case, several of the fitting parameters 
were obtained independently from fits to the numerical 
calculation thereby raising the credibility of the fits. The 
remaining parameter values (i?c2i '^i ^^d ^c) correspond 
well to those found for these materials by other means. 
When applied to another data set, |2J] Eq. (|^) was not 
found to accurately describe it. We obviously can not 
completely rule out that the results in the last reference 
are the correct ones and that the other two are wrong, but 
it seems unlikely to us. We also found that Eq. (0) could 
not describe the YBCO specific heat data of Ref. pl[ | over 
the whole temperature [80K:I00K] but that very good re- 
sults for five of the fields could be obtained if the region 
above the peak is excluded from the fit. The question of 
why this is brings us now to a broader discussion of the 
applicability of both the 3D specific and magnetization 
functions of Ref. g to YBCO data. 

As we discussed in the Introduction, several factors 
could explain possible discrepancies between the func- 



tions and the data. The most significant would be the 
failure of GL-LLL theory to describe the YBCO data. 
We believe that the evidence is against this because it 
has been shown elsewhere (as mentioned in the Intro- 
duction) that YBCO data scales according to this theory. 
Another possible reason would be discrepancies between 
the functions of Rcf. M and exact GL-LLL theory. We 
tend to discount this for two reasons. First, the 2D func- 
tion was found to have significant success in describing 
the more anisotropic HTSC materials and is known to be 
in excellent agreement with numerical simulations of the 
2D GL-LLL theory. Secondly, we note the exceptional fit 
of Eq. (P to the numerical data (Fig. |l|a) which provides 
evidence that this equation is an accurate description of 
the GL-LLL theory. Rather, we believe that the short- 
comings arise because Eqs. (^ and (|^) are 3D functions 
and it has been shown that, while the YBCO is the least 
anisotropic of the major HTSC materials, 2D signatures 
are present [^. As mentioned above, we have tested this 
by successfully fitting the 3D LLL function to the tem- 
perature range of the specific heat data which is believed 
to be 3D. Further adding to our conclusion here is prior 
evidence for 2D behavior in YBCO and LBCO through 
specific heat LLL scaling. ^,0 

The reason why the shortcomings appear to affect 
the specific heat and not the magnetization appears to 
be the following: the specific heat has a more compli- 
cated behavior than the magnetization, which is mono- 
tonic. However, the possibility that the approximations 
involved in the computation of Eqs. (^ and (||) are in- 
adequate in the region above the peak cannot be con- 
clusively ruled out at this point. To try to remedy the 
specific heat problem by splicing together the 2D and 3D 
functions in the appropriate temperature ranges to de- 
scribe the YBCO data would introduce so many fitting 
parameters that any fit would be of little value. It is 
also possible to use the quasi-2D functions of Ref. |9| but 
these also have a large number of fitting parameters be- 
sides being much less tractable than the pure 2D or 3D 
functions. 

The importance of scaling functions, as contrasted to 
mere scaling variables, has been demonstrated in this pa- 
per by showing the low information content of 3DXY 
scaling without the associated scaling functions. It was 
shown in two cases (Section |IIED how curves from an 
equation calculated using LLL assumptions could be de- 
scribed by 3DXY theory. This has consequential ramifi- 
cations for the significance of 3DXY theory and what it 
means to find that data scales according to 3DXY theory. 
In spite of the generality of 3DXY theory, it was further 
demonstrated that it does not describe the finite field 
specific heat data of the YBCO samples in the critical 
region. 
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